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ABSTRACT 

There is a recurring question in solar physics about whether or not electric currents are neutralized in ac¬ 
tive regions (ARs). This question was recently revisited using three-dimensional (3D) magnetohydrodynamic 
(MHD) numerical simulations of magnetic flux emergence into the solar atmosphere. Such simulations showed 
that flux emergence can generate a substantial net current in ARs. Another source of AR currents are photo- 
spheric horizontal flows. Our aim is to determine the conditions for the occurrence of net vs. neutralized 
currents with this second mechanism. Using 3D MHD simulations, we systematically impose line-tied, quasi¬ 
static, photospheric twisting and shearing motions to a bipolar potential magnetic field. We find that such flows: 
(1) produce both direct and return currents, (2) induce very weak compression currents — not observed in 2.5D 
— in the ambient field present in the close vicinity of the current-carrying field, and (3) can generate force-free 
magnetic fields with a net current. We demonstrate that neutralized currents are in general produced only in the 
absence of magnetic shear at the photospheric polarity inversion line — a special condition rarely observed. 
We conclude that, as magnetic flux emergence, photospheric flows can build up net currents in the solar atmo¬ 
sphere, in agreement with recent observations. These results thus provide support for eruption models based 
on pre-eruption magnetic fields possessing a net coronal current. 

Subject headings: Magnetohydrodynamics (MHD) / Sun: corona / Sun: coronal mass ejections (CMEs) / Sun: 
flares 


1. INTRODUCTION 

Current-carrying magnetic fields are an essential ingre¬ 
dient for the generation of flares and coronal mass ejec 
tions ( CMEs) in the sol ar atmosp here (e.g., Rust et al 

rjEoTuT 


1994 ; !Schrijver et ali2005klShibata & Magaral201 lt lAulanierl 


2014ft . Indeed, such non-potential magnetic fields store the 
free magnetic energy that powers these phenomena. What re¬ 
mains controversial, though, is whether and how a net electric 
current, here meant to be non-zero if integrated over one pho¬ 
tospheric magnetic polarity, is formed in the source regions of 
these phenomena. 

These questions arise from different theoretical arguments 
according to which electric currents sh ould, or should not, be 
neutralized in active regions (ARs; e.g., lMelrosell991l:|Parkerl 
11996ft . The answer to these questions may have critical con¬ 
sequences for several theoretical flare and CME models, as 
well as for pre-eruptive magnetic fields, developed fro m mag¬ 
netic configurations containing a net current (e.g., lEowll 19771: 


Tend & Kunerusll 1 978t Martenslll987UTitov & Demoulinl 


1999tlKliem & TorolJ2006 : lDemoulin & Aulanieil2010ft . In- 

deed, full neutralization implies that there js no net current in 
an AR. In these circumstances. lForbesl(l2010h pointed out that 
the eruption mechanism of these models may be inhibited. 
Their relevance may therefore be questioned if ARs currents 
are in fact neutralized. 

The question of current neutralization in ARs derives from 
the fact that the current flowing in isolated, confined mag¬ 
netic flux tubes consists of two parts: the so-c alled direct and 
return currents (lMelroselll991t fParkeri 1 199~6i) . In magneto¬ 
hydrodynamics (MHD), the direct (or main) currents refer to 
the electric currents that are expected from the chirality of a 
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twisted/sheared magnetic flux tube. For a flux rope, the direct 
currents are flowing in the central part of the twisted flux tube 
while the return currents are flowing around them (e.g., see 
Figure 3 of lMelroselll991ft . These return currents shield the 
ambient magnetic field from the direct currents. 

ARs currents are believed to be built up by two main 
mechanisms: (1) the emergence of current-carrying mag¬ 
netic flux-tubes from the solar convec t ion zone (CZ) 
into the corona (e.g.. iLeka et all Il996t jMorenoHnsertis 
1997 : iLongcope & Welschl l2000t ICheung & Isobe 

2014), and (2) the stressing of the coronal magnetic 
field by sub-photosphe ric and photospheric horizon¬ 
tal flows (e.g., [McQym pnt& E isherl 119891, 


T99TI IKlimchuk & Sturrockl 1 1 9921 ITorok & Klieml 120031 
Aulanier et al.1120lofu 

It has been argued that both mechanisms should in princi¬ 
ple produce neutralized currents. Mechanism (1) is believed 
to be associated with the rising through the CZ of confined 
magnetic flux tubes (e.g., iParkerln~955b iFanl 12009!) . For in¬ 
stance, let us consider the simplified case of a twisted flux 
tube carrying an electric current, I, in cylindrical geometry. 
The confinement of the flux tube to a finite cross section of ra¬ 
dius, R, requires that the total current it carries must vanish for 
r > R (see Appendix lAli. Such a twisted flux tube possesses 
a non-zero internal current-density, j. Therefore, the flux tube 
must contain a second type of currents which neutralizes the 
core/direct currents: i.e., having the same total strength but 
flowing in the opposite direction (and often assumed to be a 
surface current). Based on the simplified assumption of full 
emergence of confined magnetic flux tubes, one may then ex¬ 
pect that mechanism (1) would transport neutralized currents 
into the solar corona, thus generating a current-neutralized 
AR. As for mechanism (2), localized (sub)-photospheric hor- 
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izontal flows transfer twist and shear to the magnetic field in a 
finite coronal volume. This field will typically inflate, induc¬ 
ing currents through compression also in the ambient mag¬ 
netic field. Since the (sub)-photospheric driving volume is 
finite, one may expect that the changes in the coronal field 
also remain restricted to a finite volume. If true, a complete 
shielding of the generated currents, i.e., neutralized currents, 
would be implied. 

Observationally, the normal/vertical component of the elec¬ 
tric current density, j z , can be derived by applying Ampere’s 
law to photospheric vector-magnetograms (see Equation (0}). 
Despite the various uncertainties and difficulties, the mea¬ 
surements of photospheric transverse magnetic fields are be- 
coming more and more reliable (e.g. , iL eka et al. 1119961: 
Metcalf et all l2006t IWiegelmann et all l2006t IGosain et all 
20141) . For this reason, increasing attention has been paid to 
deriving the properties of currents in solarARs^ and to test¬ 
ing their degree of neutrali z ation (e.g.. iWilkin son et all 1992; 
Leka & Skumanichl 119991 IVenkatakrishnan & Tiwaril 2009 : 
Sun et al.ll2012l) . Some recent observational studies report the 
presence of direct a nd return currents in each magnetic po¬ 
larity of ARs (e.g.. lWhea tlandlf2 000t iRavindra et al.l 120111 
iGeorgoulis et al.l2012tIGosain et al.l2014l) . These studies find 
both types of ARs, i.e., some with neutralized currents, and 
some with a net current. 

The indications for the existence of ARs with a net current 
are at variance with theoretical arguments invoked in favor 
of current neutralization. Considering their past and present 
limitations, the relevance of the observational measurements 
has thus bee n questioned (e.g..lPark er 1996; see also the In¬ 
troduction of lGeorgoulis et al .1120 1 2l) . On the other hand, the 
arguments in favor of current neutralization may well be over¬ 
simplified. For instance, they usually do not consider possible 
effects that become relevant in a fully three-dimensional (3D) 
geometry or in the case of partial magnetic flux emergence. 
This has thus led to a long-lasting debate about whether or 
not a net current c an exist in ARs (see e.g., lMelrosel|l9911 
[T99llParkeil[T996h . 

Numerical MHD simulations provide a useful alterna¬ 
tive for addressing this problem. However, the neutraliza¬ 
tion of electric currents has barely been analyzed with this 
tool. A few MHD simulations reported the presence of 
both direct and return currents generated by photospheric 
line-ti ed motions applied to initially potential co ronal fields 
(e.g.. lAu lanier et ~al)120051; iDelannle et al.l[20081) . Yet, only 
ITorok & Klieml (120031) quantified the associated degree of 
current neutralization. For the case of a twisted flux tube, 
they found that the neutralization only occurs when the pho¬ 
tospheric motions do not extend to the polarity inversion line 
(PIL), so that no magnetic shear is built up at the initially un¬ 
sheared PIL. 

ITorok et all (120141) were the first to revisit the question of 
current neutralization by means of 3D MHD simulations of 
magnetic flux emergence. In their experiment, the authors 
modeled the emergence of a buoyant magnetic flux rope into 
a stratified, plane-parallel atm osphere in hydrostatic equi¬ 
librium (see iLeake et al.l [20131) . For that purpose, a sub- 
photospheric magnetic flux rope containing neutralized cur¬ 
rents was considered. It was found that a complex redistri¬ 
bution of the initially sub-photospheric direct and return cur¬ 
rents occurs in the vicinity of the photosphere. This subtle 
redistribution led mainly to the emergence of the initial direct 
currents (see Figure 3b of lTorok et al.ll20141) . causing the de¬ 
velopment of a strong net current in the corona. This net cur¬ 


rent was associated with the development of a strong magnetic 
shear along the PIL, and s ome non-force-free return currents 
(see Figure 5 of lTorok et al.l2014h . These results indicate that 
the emergence of a current-neutralized magnetic flux tube can 
lead to the generation of a n et current in ARs, as suggested by 
iLongcope & Welsriil (120001) . 

In the present study, we pursue the work initiated by 
ITorok etaTI (120141) . analyzing here the distribution and neu¬ 
tralization of currents generated by photospheric horizontal 
flows. We perform a parametric analysis of 3D, zero-/?, MHD 
simulations of photospheric twisting and shearing motions 
imposed on a bipolar potential field. Section [2] describes the 
main set-up of our numerical models. The results of our para¬ 
metric study are presented in Section 0] for the twisting mo¬ 
tions and in Section 0] for the shearing ones. A discussion 
and interpretation of our results is provided in Section 0 Our 
conclusions are summarized in Section[6] 

2. NUMERICAL MODEL 

2.1. Equations, numerical scheme, and boundary conditions 

The numerical simulations described in this paper were 
performed in cartesian geometry, using the Observationally- 
driven High-o rder scheme Magnetohydrodynamic code 
(OHM; see lAulanier et al.ll2005l) . We use the code in its zero- 
fj version in which the mass density, p, the fluid velocity, u, 
and the magnetic field, B, are advanced in time according to 

dp 

^ = -V • (pu), (1) 

— = - (u • V) u + (V x B) x B /(pop) + vAu, (2) 

ot 

O T> 

— =Vx(uxB) + /7AB, (3) 

at 

where v and q are diffusion coefficients, v is a pseudo¬ 
viscosity, and // is the electrical resistivity. These diffusions 
coefficients are used to limit the development of sharp dis¬ 
continuities that may develop at the scale of the mesh and lead 
to quickly-growing nume rical instabilities (see lAulanier et al] 
120051; Uanvier et al .120 131) . The solenoidal condition (V B = 0; 
a discussion on its very weak magnitude is provided in Ap¬ 
pendix [B]) and the current density are not calculated in the 
code. The latter is derived from Ampere’s law 

j = —V x B . (4) 

Eo 

Equations £0-0} are solved in non-dimensionalized units, 
using /To = 1. The velocities are expressed in units of the 
initially uniform Alfven speed, c A (t = 0) = 1. The time 
unit is given by the Alfven time, t A — 1, which corresponds 
to the travel time of an Alfven wave over a distance d = 1. 
The diffusion coefficients are prescribed in terms of uniform 
characteristic speeds, Uy and u n , such that v = Uy/lijy and 
q - u rl lijy, where lijy is the smallest grid-spacing at point 
(xpypik)- For all simulations, we set Uy = 15 %c A (t - 0) 
and Ujj = 1 ,5%c A (t = 0) in the coronal volume (the diffusion 
parameters are set to zero at the photospheric line-tied bound¬ 
ary). 

All simulations are performed in the same domain cover¬ 
ing x xy x z ~ [-9.1,9.1] 2 X [0,30], discretized on an non- 
structured mesh of n x x n y x, n 7 = 231 3 points. In the sub¬ 
domain .r x v = [-1.5,1.5] 2 , the mesh is set uniform along the 
x and y directions, with mesh intervals d™ n — of™ n = 0.02. 
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This choice allows us to employ the higher mesh-resolution 
in the region where the strongest gradients of electric cur¬ 
rents are generated. Outside of this region, the mesh is non- 
uniform in both the x and y directions, with mesh intervals 
defined such that d'+ l /d' x — d {' 1 /d{ = 1.091. The mesh is 
non-uniform all along the z direction, with mesh intervals de¬ 
fined by d k z +1 /d* = 1.013 and df n = 0.02. With these set¬ 
tings, the largest mesh intervals are c/ lnax = 0.41 along z and 
= 0.65 along both the x and y directions. 

We use the same initial conditions for each numerical sim¬ 
ulation with an initially potential magnetic field constructed 
by placing two fictitious opposite magnetic charges below the 
photosphere. The positive and negative charges are placed at 
r ± = (0, +1, -1) respectively. The resulting potential field is 


r - r + r - r 

B />( r ) = go |„ _ |3 ~ qo 


r - r_| 3 ’ 


(5) 


where qo(> 0) is the strength of the two magnetic charges. For 
all simulations, we use qo — 1. The corresponding magnetic 
field is displayed in Figure |T| The initial plasma density is set 
to po = B 3 /(poc^), where ca is the initially uniform Alfven 
speed. 

The top and all side boundaries of the numerical domain are 
open (further details^ on th ese boundar y conditions are given 
in Section 2.5 of lAulanier et al.l 12005 ). Line-tied boundary 
conditions are prescribed at the photospheric, or z = 0, plane 
to build-up twisted and sheared magnetic fields. 


2.2. Photospheric twisting motions 

The first type of photospheric driving is applied to twist 
the initial potential magnetic configuration along the isocon¬ 
tours of the photospheric vertical magnetic field, B z (z = 0). 
This modifies the transverse components of the initial poten¬ 
tial field while preserving its photospheric flux distribution, 
B-(z = 0 ,1 - 0). The expression of the twisting velocity field 
is 


dip 

K x(z = 0 ) = Mo — , 
ay 

(6) 

dip 

u v (z = 0 ) = -m 0 — , 

OX 

(7) 

ii 

o 

II 

o 

(8) 


where uq is a free parameter that controls the maximum speed 
of the driving, an d ip is a normaliz e d, time-dependent po¬ 
tential (see e.g., | Amari et all 1 19961 iTorok & Klieml 120031 
lAulanier et al] 120051) . This potential depends on B z (z = 0) 
such that 

_ B\{z = 0) ( B 2 z (z = 0)- 2T"- 2 fe = 0) ) 

* = o) exp [ tiB'r Hz = o) j ’ 

This twisting profile generates two vortices centered on 
±B" vdx (z. = 0) respectively. The size of the vortices is con¬ 
trolled by the free-parameter it*. 

All twisting profiles were applied with a maximum driving 
speed uq = ().()2c A (t = 0). Figure |T] shows the profiles of 
ujx = 0,y) for all four cases considered in our study, and 
referred to as T{ 1; 2; 3; 4). The corresponding parameters are 
listed in Table[t| The left panels of Figure [2] display the 3D 
distribution of the magnetic field for T3. 

While the twisting boundary driving analytically preserves 
the initial distribution of the photospheric vertical magnetic 
field, small numerical errors eventually deform it on the long 
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Twisting profiles 




Fig. 1. — Top: initial magnetic field configuration of the shearing and twist¬ 
ing simulations. The synthetic magnetogram is represented at the z = 0 
plane with superposed ±[0.25; 0.50; 0.75] isocontours of B z (solid purple and 
dashed cyan) and selected field lines (red). The green line shows the polarity 
inversion line. Center: photospheric (z = 0 ) B r and u x profiles in the y di¬ 
rection for the twisting motio ns (s ee Section ET2l . Bottom: same but for the 
shearing motions (see Section [23). 
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4 cases 4 cases 2 cases 

Fig. 2.— Top (top row) and 3D (bo ttom row) views of the magnetic field g ener ated by t he photospheric line-tied motions. Left: twist case T3 (Equations J6]- 
[9}). Center: shear case S2 (Equation ) 101 ). Right: shear case SP2 (Equations i 11 and )12t ). Current-carrying (potential) magnetic field lines are in color (black). 


TABLE 1 

Values of the parameters for the velocity profiles 


Run 

itw 

u 

£shl 

£shl 

yo 

yi 

n 

T1 

0.1 

- 

- 

- 

1 

- 

_ 

T2 

0.5 

- 

- 

- 

1 

- 

- 

T3 

1 

- 

- 

- 

1 

- 

- 

T4 

5 

- 

- 

- 

1 

- 

- 

SI 

- 

0.17 

- 

- 

1 

- 

- 

S2 

- 

0.32 

- 

- 

1 

- 

- 

S3 

- 

0.55 

- 

- 

1 

- 

- 

S4 

- 

0.77 

- 

- 

1 

- 

- 

SP1 

- 

- 

0.24 

0.24 

- 

0.5 

0.5 

SP2 

- 

- 

0.24 

0.63 

- 

0.5 

0.3 


Note. — All values are in non-dimensional units 
(normalized by half the distance between the two pho¬ 
tospheric magnetic polarities). 


run. To ensure that B- is preserved in time for the twisting 
simulations, we numerically reset B.(z = 0, t) to B~(z = 0, t = 
0) at each time step in the course of the runs. 

2.3. Photospheric shearing motions 

We consider two types of photospheric shearing motions 
along the x direction, u(z = 0) = (u x , 0,0). The first is 
centered on the strongest magnetic field of both photospheric 
magnetic polarities. The second is confined to the weak field 
surrounding the PIL to primarily shear it (e.g., similarly to 
[Antiochos et al.lfl994l) . 

The first type of shearing motions (curves S{ 1; 2; 3; 4} in 


Figure H} is given by 




( (y-yo) 2 l 


(y + yo) 2 

MQ AO 

exp 

r2 

^ ^ sh 2 

- exp 



( 10 ) 

where controls the width of the shearing, yo corresponds to 
the y-position of the maximum velocity (shearing center), and 
a o is a constant used for normalization of the Gaussian-profile. 
The resulting shearing is invariant along the x direction. 

The second type of shearing motions (curves SP{ 1; 2} in 
Figure Q} is given by 


u x (z = 0) = M 0 fli 


exp 


A 


- exp 


for |y| < yi, and 


(y + yi) 2 Y 

a 

G/ii 

( 11 ) 


ll x (z = 0 ) = Mq «2 


exp 


(y - V2) 2 ' 

A 
b shl 


- exp 


(>’ + V2) 2 ' 

A 

±sh2 / 


( 12 ) 

for |y| > y 1 . Choosing £^1 + £ s h 2 allows one to have a broader 
sheared region either for |y| < y 1 or |y| > y 1 . The parameter y 2 
is computed to ensure the continuity of the shearing profile at 
lyl = yu where u x (\y\ =yuz = 0) = u 0 . 

In this paper, we consider four shearing profiles for the first 
model, referred to as S{ 1; 2; 3; 4}, and two shearing profiles 
for the second model, referred to as SP1 and SP2. The corre¬ 
sponding parameters are listed in Table [Q 
As for the twisting simulations, all shearing profiles were 
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applied with a maximum driving speed uq - 0.02 cfit = 0) to 
ensure a quasi-force-free evolution of the magnetic field. The 
applied profiles are displayed in FigureQ] The middle and left 
panels of Figure[2]display the 3D distribution of the magnetic 
field for S2 and SP2, respectively. 


2.4. Ramp function 

We apply all photospheric line-tied motions using a ramp 
function to smoothly bring the system from rest to a constant- 
velocity photospheric driving, such that 


u (x,y,z = 0 ,t) = y(f) u (x,y,z = 0), (13) 

where u(x, y, z = 0) is defined in Sections 12.21 and 12.31 The 
ramp function, y(t), is given by 


jit) = - tanh 


2(t - t m ) 



(14) 


where t m corresponds to the time at which the middle of the 
ramp function is reached, and th w corresponds to the half¬ 
width of the ramp time. In our runs, we fix t m = 1 5t A and 
thw — 5 tA- With these settings, the photospheric driving starts 
after » 10 Alfven times, allowing the system to reach a good 
numerical equilibrium before the acceleration begins. The 
constant-velocity photospheric driving is reached after « 20 
Alfven times. 


3. PHOTOSPHERIC CURRENTS INDUCED BY TWISTING MOTIONS 

In cylindrical geometry, any twisting motion based on a 
twist function that falls to zero at a finite radial distance should 
generate a twisted flux tube formed of a core of direct cur¬ 
rents, fully surrounded by a shell of return currents that ex¬ 
actly neutralize the direct currents (see Appendix lAli. 

The transposition of these results to 3D geometry has of¬ 
ten been used to argue that compact photospheric twisting 
motions shoul d lead to the generation of fully neutralized 
currents ( e.g ., [Melrosel 1 199 It iParkerl11996b . However, such 
a transposition disregards the fact that a 3D magnetic flux 
tube may not systematically have a cylindrical analogue. It 
is therefore not obvious that the properties of electric cur¬ 
rents expected from a simplified, cylindrical geometry may 
still hold in a more complex 3D one. 

To address this issue, we analyze the electric currents gen¬ 
erated by compact, photospheric twisting motions (see Sec¬ 
tion [2T2j>. 


3.1. Photospheric distribution of vertical currents 

Figure [3] displays the photospheric vertical current den¬ 
sity, j z , for the twisting runs, a few Alfven times before 
each numerical simulation terminated (due to the develop¬ 
ment of a numerical instability caused by very sharp gradi¬ 
ents). Each twisting run induces the generation of a core 
of negative-direct currents surrounded by a shell of positive- 
return currents, just as in cy lindrical geometry. As reported by 
iTorok & Klieml (120031) and iDelannee et al.l (120081) for similar 
models, we find that the direct currents are overall stronger 
and more compact than the return currents. In addition, we 
notice that the distribution of currents always exhibits the 
same type of strong azimuthal asymmetry (except for Tl) that 
does not occur in cylindrical geometry. 

This gener ic _ propert y was al r eady explained by 

lAulanier et al.l (120051) and lTitov et al.l (120081) . and is typical 
of the 3D loop-geometries analyzed in this paper. This 
asymmetry is an effect of field-line length resulting from 
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Fig. 3.— Photospheric (z = 0) current maps, j z , at the end of each simu- 
lati ons f or the symmetric twisting of two opposite magnetic polarities (Sec- 
tion l2.2l , for T{ 1; 2; 3; 4}. White and black display positive and negative cur¬ 
rents respectively. Values are saturated at ±1.5. Solid purple and dashed cyan 
lines represent ±0.25 isocontours of B z . 


the flux tube curvature. In cylindrical geometry (r, 6 , z), the 
equations of magnetic field lines and electric current density 
lead to 


rO 

B e =—B z , 


. 1 drBg Bg 

Jz= -oc — cx 

r or r 


1 ’ 


05) 

(16) 


where ® and l are the field-line twist and length, respectively. 
For a given B z isocontour <t> is fixed. It then follows that 
shorter field lines (i.e., smaller /) possess a stronger electric 
current density (Equation ([16])). For our 3D curved flux tube, 
the same amount of twist, <t>, is transferred to each field line 
of any given isocontour of B z . The previous considerations 
thus imply that stronger currents develop at the footpoints of 
the shortest field lines of any B z isocontour. This creates an 
asymmetry that is amplified by the faster expansion of the 
larger field lines. The advection of this asymmetry by the 
photospheric motions is responsible for the swirling pattern 
exhibited by the electric current distribution. This effect is 
extremely weak for Tl because the twisting vortices are so 
narrow that the twisted field lines have a very similar length. 

Among several differences between each simulation, one is 
the extension of the currents close to the PIL. In particular, for 
the cases T{ 1; 2}, the twisting vortices are so narrow that the 
distribution of current is mainly localized well inside the iso¬ 
contour \B Z \ = 0.25. On the contrary, for T{3;4), the vortices 
are so broad that the distribution of current extends to the PIL. 


3.2. Evolution of the total direct and return currents 

We now analyze the curves of integrated currents by com¬ 
puting 


'? = 


-f k 

Jy> 0 


dS, 


(17) 
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Fig. 4. — Evolution of the electric current, I z , at z = 0 in the positive po¬ 
larity for the photospheric twisting motions. Top: negative-direct (7? irect ) and 
positive-return f/' clllr ") currents (respectively displayed in solid and dashed 
lines). Bottom: neutralization ratio. 


where X refers to the total direct or return current within 
the positive magnetic polarity (i.e., y > 0). Because the 
direct/return currents are negative/positive, we compute the 
total direct/return current by extracting the negative/positive 
current density in the positive magnetic polarity. 

The temporal evolution of the integrated direct and return 
currents is presented in the top panel of Figure 0] Due to its 
small vortex size, the T1 case shows the development of only 
weak direct and return currents. For the three other cases, the 
vortices are broader and stronger direct and return currents 
develop. 

Figure 0] shows that the absolute value of the integrated di¬ 
rect current exhibit a monotonic rise (regardless of the vortex 
width) and the strength of the direct current increases with 
the width of the twisting vortex. The integrated return current 
also presents a monotonic behavior. However, the strength of 
the return current only increases with the vortex width up to 
(tw < 0-5. 

Increasing the vortex width builds up twist in a larger vol¬ 
ume and generates a higher direct current. The return current 
is affected in a more complex manner, because the boundary 
between the direct and return currents is pushed closer to the 
PIL. Away from the PIL, increasing the vortex also implies 
that longer magnetic field lines are being twisted. As dis¬ 
cussed by lAulanier et al. ( 2005 ). the resulting fast expansion 
of these field lines limits their current density. 

We therefore conjecture that the increase of total electric 
current (|/? Ilect | + |/j etum |) with the vortex size, is caused by 
a complex competition of three mechanisms: (1) an increase 
due to the transfer of magnetic twist in a broader region, (2) a 
saturation/decrease caused by the approach of the line of cur¬ 


rent reversal to the PIL, and (3) a saturation/decrease induced 
by the fast expansion of the magnetic field lines. 


3.3. Evolution of the neutralization ratio 

In order to quantify the neutralization of electric currents 
generated by photospheric line-tied motions, we define the 
neutralization ratio, J neu t., as 


I 


neut. 


yreturn 


/direct 

1 z 


(18) 


where 7* rect and / rulurn are computed from Equation ( |T7i The 
neutralization ratio is 1 for fully neutralized currents and 0 for 
a magnetic field solely containing direct currents. 

The bottom panel of Figure0]shows the temporal evolution 
of the neutralization ratio for the twisting runs. Our goal is 
to analyze the neutralization of currents for different spatial 
profiles of photospheric line-tied motions. In the following, 
we therefore restrict our analysis to the period during which 
the system evolves in response to a stationary boundary driv¬ 
ing, i.e., for 1 > 20t A . A brief discussion of the results for 
the transition phase, 1 e [0; 20] Lt, is nonetheless provided in 
AnnendixICl 

We find that the electric currents remain fully neutralized 
(J neut . = 1) during the run with the narrowest vortices (Tl), 
and nearly neutralized (I~ n eut. ~ 0.98) for the T2 run. On the 
contrary, the two runs with the larger vortices, T{3; 4}, exhibit 
a strong departure from neutra lizatio n. This confirms the ear¬ 
lier results of lTorok & Kliemi (120031 ) that were obtained with 
different numerical settings. 

Figure 0] further indicates that the constant boundary driv¬ 
ing phase is associated with a slow and weak decrease of the 
neutralization ratio for run T3. On the contrary, the neutral¬ 
ization ratio of run T4 presents a weak increase followed by 
a saturation. Such behaviors result from a non-trivial com¬ 
bination of the increase of currents with magnetic twist and 
the saturation of currents caused by the fa st ex pansion of the 
largest field lines, as discussed in Section 13.21 Note that for 
T4, two additional effects are likely involved in the evolution 
of the neutralization ratio: (1) the twisted flux tube starts to 
leave the numerical domain (as suggested by the opening of 
some field lines, not shown here), and (2) the flux tube enters 
a super-exponential growth phase (cf. the equi librium curves 
in iTorok & Kliemi 12003b lAulanier et al.l l2005h . leading to a 
fast expansion of the more twisted field lines compared with 
the less twisted ones, i.e., for the direct currents. This lim¬ 
its the growth of the direct current more efficiently than that 
of the return current, which could explain the observed weak 
increase of the neutralization ratio. 

In contrast to cylindrical symmetry (AnnendixlAli. Figure0] 
shows that twist profiles do not systematically generate neu¬ 
tralized currents in fully 3D twisted flux tubes. Hence, the 
condition for current neutralization derived in 2.5D geome¬ 
try cannot be directly transposed to a fully 3D loop geometry 
such as the one considered in this paper. Indeed, we will show 
in Section[5]that (1) the condition for current neutralization in 
3D is more subtle than in 2.5D, and (2) net currents are asso¬ 
ciated with 3D coronal fields that do not have any analogue in 
2.5D. 


4. PHOTOSPHERIC CURRENTS INDUCED BY SHEARING MOTIONS 

In this section, we study the electric currents gen erate d by 
compact photospheric shearing motions (see Section 03}. 
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Fig. 5.— Photospheric current maps as in Figur e [T| f or the symmetric shear¬ 
ing of two opposite magnetic polarities (Section 123F Top/middle: shearing 
centered on B™ ax for S{ 1; 2; 3; 4}. Bottom: shearing centered on the weak 
field surrounding the PIL for cases SP{ 1; 2}. The color coding is the same as 
in Figure[3|with saturation at ±0.5. 

4.1. Photospheric distribution of vertical currents 

Figure [5] displays the photospheric vertical current density, 
j z , for the shearing runs. We find that each run presents both 
negative-direct and positive-return, force-free currents. Their 
sign are identified from the sign of magnetic helicity trans¬ 
ferred to the system. In particular, the chosen shearing mo¬ 
tions produce negative magnetic helicity, and hence, negative- 
direct currents. 

For almost all simulations, the direct and return currents 
have a similar spatial distribution with similar intensities. Two 
cases, however, do not display the same pattern: S3 and S4 
mostly possess direct currents. For these two runs, the return 
currents are much weaker than the direct ones ( e.g ., « 7 times 
weaker for S3). 

The return currents are, in general, expected for localized 
shearing motions. This is qualitatively described in Figure [6] 
with the right-hand rule relating the local magnetic shear with 
the direction of the electric current. Localized shearing mo¬ 
tions generate a curl of the magnetic field that changes sign 
at the line of stron gest magneti c shear. However, contrary to 
the drawing of iMelrosel ( 1991 ). the line of current reversal 
may not systematically correspond to that of the maximum 
velocity. For instance, the maximum photospheric velocity 
occurs at |y| = yo = 1 for the S{ 1; 2; 3; 4} runs. Nevertheless, 
only S{ 1; 2} possess a current reversal at |v| ~ yo at the pho- 



-2 X 2 -2 X 2 


Fig. 6. — Examples of sheared magnetic field lines and current direction 
generated by photospheric shearing motions for the case S2. The thick blue 
(respectively green) field line is anchored within the direct (respectively re¬ 
turn) current. The thick black field line is anchored at the line of current 
reversal (thin yellow lines). The orange arrows and hands show the direction 
of current density, j, on both sides of the line of current reversal, as inferred 
from the right-hand rule. The arrows on the field lines indicate the magnetic 
field direction. The color coding and saturation are the same as in Figure[5] 

tosphere, while it occurs at |y| < 0.5 for S3 and at the PIL for 
S4. This is because the current distribution depends not only 
on the shearing velocity profile but also on the 3D shape of 
magnetic field lines. 

As for the twisting runs, we find that one significant dif¬ 
ference between each shearing simulation is the extension of 
the currents close to the PIL. In particular, the shearing region 
of SI and S2 is so narrow that their distribution of currents 
is strongly localized (in the v-direction) within the isocontour 
\B Z \ = 0.25. On the contrary, the shearing region is so broad 
for S3 and S4 that the distribution of currents extends to the 
PIL. Finally, both the shearing profile and electric currents 
extend to the PIL for SP1 and SP2. 

4.2. Evolution of the total direct and return currents 

The top panel of Figure [^presents the temporal evolution of 
the integrated direct and return currents in the positive mag¬ 
netic polarity. All direct current curves show a similar mono¬ 
tonic rise in absolute value with a higher intensity for broader 
shearing regions. By contrast, the evolution of the return cur¬ 
rents varies with the width of the shearing region. For S{ 1; 2}, 
the return current monotonically increases with time. For S3, 
the intensity of the return current increases and then smoothly 
decreases, while for SP{ 1; 2} the simulation terminated too 
early to draw conclusions. The S4 run does not show any 
return current beca use the line of current reversal occurs at 
the PIL ( cf. Section [4~il . 

As shown in Appendix[D] the applied shearing motions can 
generate two evolutionary phases for the current density of 
magnetic field lines: (1) a first phase of increase with mag¬ 
netic shear, and (2) a phase of decrease caused by a fast elon¬ 
gation of the field lines. For similar shearing velocities, the 
shortest field lines should be the first to experience (2), be¬ 
cause this effect is more pronounced for field lines that align 
more rapidly with the PIL (as explained in Appendix iDli. For 
all shearing simulations, the return currents are associated 
with the shortest magnetic field lines. We thus argue that it 
is the late elongation of these field lines that is responsible for 
the decrease of return current observed for S3. 

Finally, we note that increasing f s ya for the second shearing 
profile (i.e., SP2), allows one to modify the shear profile of the 
field lines associated with the direct current. In particular, the 
direct current is distributed over a broader region. However, 
the amount of direct and return current is preserved. 

















8 


Dalmasse et al. 



Fig. 7.— Evolution of electric current as in Figure [4] for the photospheric 
shearing motions. 


4.3. Evolution of the neutralization ratio 

The bottom panel of Figure [7] displays the temporal evolu¬ 
tion of the neutralization rat io fo r each shearing run. For the 
reasons explained in Section l3T3l we focus our analysis on the 
neutralization during the constant photospheric driving phase 
(, i.e ., t > 20 ? a )- 

The currents remain fully neutralized (J" neut = 1) for SI and 
S2. On the contrary, the two runs with the broadest shearing 
widths (S{3; 4}) exhibit a strong departure from neutralization. 
For S4, the neutralization ratio vanishes as a consequence of 
the absence of return current in the simulation (becaus e th e 
line of current reversal occurs at the PIL; cf. Section 14.11 1. 
For SP{1;2}, the shearing motions applied close to the PIL 
also generate a strong net current. Both runs lead to the same 
degree of current neutralization. 

We notice that the curves of neutralization ratio of the S3 
and SP{ 1; 2} runs all present a comparable evolution. They all 
show a continuous decrease. Such a decrease indicates that 
the rate of current build-up is stronger for the direct current 
than for the return current. As previously mentioned (Sec¬ 
tion B3 >, the effect of field line elongation more strongly af¬ 
fects the shortest field lines. One then expects a lower rate of 
current build-up for the return current than for the direct one. 
This could then explain the continuous decrease of neutraliza¬ 
tion ratio observed for the SP{ 1; 2} and S3 runs. 

Finally, the decreasing behavior of the neutralization ratio 
of S3 and SP{ 1; 2} indicates that, in a system driven by station¬ 
ary shearing motions, the neutralization state of the system is 
not solely determined by the spatial properties of the shear¬ 
ing motions, but also depends on the amount of accumulated 
shear and inflation. 

5. DISCUSSION 



-3 -2 -1 X I 2 3 -3 -2 -1 X 1 2 3 


Fig. 8. — Saturated photospheric j- (gray shading) showing the currents as¬ 
sociated with the compression of the ambient held caused by the inflation of 
the twisted (left) and sheared (right) magnetic fields. The saturation value is 
2.4 X 1CT 3 for T3 and 5 X 1CT 3 for SI (compared with Figures[2and[5]where 
the saturation value is two orders of magnitude higher). The solid-purple and 
cyan-dashed lines are |S-| = 0.25-isocontours. 

In this section, we first examine the development of weak 
compression currents in the ambient field of our simulations. 
We then discuss our results in the framework of the sheared- 
PIL/net current relation ship and confront them to the conjec¬ 
tures of iMelrose] (1199 ll) and [Parker] (119961) . 

5.1. Compression currents in the ambient field 

When analyzing the distribution of currents of our 3D 
current-carrying fields at high saturation levels, we find that 
some very weak currents develop in the ambient field rooted 
in areas where the photospheric flows vanish (as pointed by 
the red arrows in Figure[8]). These currents are typically 2-3 
orders of magnitude smaller than the currents directly gen¬ 
erated by the photospheric flows. They appear in the very 
close vicinity of the twisted/sheared fields and very rapidly 
decrease away from them. We find that these currents are in¬ 
duced by the local compression of the ambient field caused 
by the inflation of the twisted/sheared field in response to the 
photospheric flows. Such a compression of the ambient field 
cannot be reproduced in 2.5D geometry because the imposed 
symmetry forces that field to stay potential. 

We also note that the compression currents do not form a 
shell of a single sign around the return current for the twist¬ 
ing cases ( cf. T3 in Figure [8J. On the contrary, we find two 
regions of enhanced currents oriented in two specific direc¬ 
tions (as indicated by the dashed red lines in Figure [8}. This 
is caused by the development of a kink of the flux tube axis. 
The direction of the L-dashed red line corresponds to the ori¬ 
entation of the kink of the flux tube axis (as shown in Figured 
top left and more precisely by the yellow field line in Figure 
5 of lAulanier et al.ll2005l) . The kink of the axis causes a pref¬ 
erential compression of the ambient field in its direction, and 
induces a magnetic depression in the direction of the /^-dashed 
red line. 

5.2. The sheared-PIL/net current relationship 

The numerical experiments of Sections [3] and [4] show that 
photospheric line-tied motions can generate 3D force-free 
magnetic fields with different amount of current neutraliza¬ 
tion. We further analyzed and compared the simulations pre¬ 
sented in this paper to identify the origin of these various de¬ 
grees of neutralization. The results are summarized in Fig¬ 
ure [9] that displays the photospheric transverse magnetic field 
at the PIL for current-neutralized and net current cases. We 
find that all current-neutralized magnetic fields possess a PIL 
fully embedded in a potential magnetic field (Figure [9] left 
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Fig. 9. — Sheared-PIL/net current relationship illustrated for two current- 
neutralized (T1 and SI; left column) and two net current (T4 and S4; right 
column) simulations. The photospheric transverse magnetic field at the PIL, 
B r (PIL) (green arrows), is plotted over the photospheric B z (gray shading). 
The solid-purple and cyan-dashed lines are \B Z \ = 0.25-isocontours. 

column). On the contrary, we observe that (1) a net current 
develops simultaneously with magnetic shear at the PIL, and 
(2) a stronger net current is induced for simulations generat¬ 
ing a stronger magnetic shear jit th e PIL (see Fig ure [9] right 
column). In agreement with lTorok & Klieml (12003 ). and ex¬ 
tending their results to pure shearing motions, we thus find 
that a net current solely develops when magnetic shear is built 
up at the PIL. These results are also consistent with obser¬ 
vation al studies ( e.g ., iRavindra et al.ll20llt iGeorgoulis et al. 
|2012h and the investigation of flux emergence in lTorok et al. 
(120141) . In the following, we address this relationship in a 
general form. 

We consider a bipolar potential magnetic field such as the 
initial field used in our numerical simulations (Figure [TO]). We 
then consider a photospheric velocity field, u, that builds up 
electric currents at z — 0 inside of a surface, S u , of the posi¬ 
tive magnetic polarity 1 ; u vanishes outside of S u . Finally, we 
consider a closed curve, C, such that (i) C includes the PIL of 
the bipolar AR, and (ii) the surface, S , enclosed by C is much 
larger than the surface S u where electric currents are gener¬ 
ated by the photospheric motions. These two choices ensure 
that all the currents transferred to the magnetic field are fully 
enclosed by C. 

The applied photospheric velocity field generates a current- 
carrying component. By, in the bipolar magnetic field. The 
total magnetic field, B, can be decomposed as the sum of 
its pote ntial, B„, and current-carrying. By, components (e.g., 
IValori et ai1l2013l) , such that 

B = B /) + By, (19) 

1 The corresponding definition can also be done in the negative magnetic 
polarity, but this is not needed because of V ■ j = 0 and current is transferred 
to closed magnetic field lines. We therefore keep our analysis simpler by 
focusing on the positive magnetic polarity, as we did with the analysis of our 
numerical simulations. 



-5 0 5 

Fig. 10.— Schematic used for the mathematical demonstration of the 
sheared-PIL/net current relationship. The map shows a photospheric magne¬ 
togram of the initial potential magnetic field used in our line-t ied s imulations. 
The surfaces S/S u and the contour C are defined in Section PT2l The thick, 
black/white, dashed line highlights the part of C that corresponds to the PIL 
(green solid line). The solid purple, and dashed cyan, lines are B z = +0.25- 
isocontours of the magnetic field. 

V x B = V x B ; = yuoj . (20) 

In fact, B j can be further decomposed into two components 
By = Bf°+ B“ =0 . (21) 

B“ ,: ° corresponds to the main current-carrying field and it has 
non-zero values only in the magnetic field connected to S u - 
B“ =0 is a current-carrying field generated by the compression 
of the potential field resulting from the inflation of the main 
current-carrying field. This component is non-zero only in 
the ambient magnetic field, i.e., the magnetic field not con¬ 
nected to S u . Its strength is on average two orders of mag¬ 
nitude smaller than the strength of the main current-carrying 
field, B“^°, and decreases very rapidly away from the latter 
(cf Section [57Tb . For this reason and because the contour C 
is chosen such that S is much larger than S u , we can neglect 
B“ =0 in the following. 

Applying Ampere’s law, it then follows that the total elec¬ 
tric current enclosed by C is 

I=-£bj- dl 
Mo Jc 

~—(D B u#0 • dl, (22) 

Mo Jc 1 

where dl is the line element along C. Equation (l22l > can 
be further decomposed as the sum of two contributions: one 
along the part of C corresponding to the PIL, and a second 
corresponding to the remaining of C (i.e., C - PIL). This 
leads 

/ « — ( f B 11,60 • dl + f B 11,60 • dl) . (23) 

M0 \JPIL 1 JC-PIL 1 I 

The intersection of C with S u — when it exists — is limited 
to the PIL, and B“ /0 vanishes outside S u . It then follows that 
the second term of Equation (l23l) vanishes. At any fixed time. 
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the total current enclosed by C is therefore 

/ ~ — f Bf 0 • dl . (24) 

Po JpiL 1 

When the PIL is embedded within a potential magnetic field, 
i.e., S ,i is never in contact with the PIL, then B" /() vanishes 
all along the PIL. The total electric current, 7, thus vanishes 
as well, which implies that the currents are neutralized. On 
the contrary, B“*° does not vanish for a magnetically-sheared 
PIL. The current-carrying magnetic field will then contain a 
net electric current, 7 4 0 (unless the integrand in Equa¬ 
tion (124b changes sign such that the oppositely directed con¬ 
tributions cancel exactly, which could happen for the very rare 
cases of ARs containing opposite twist/shear). We thus con¬ 
clude \halforce-free net currents are inevitably related to mag¬ 
netically sheared PILs. 

The above derivation is supported by the results of the SP2 
run as compared with SP1. Indeed, both possess the same 
shearing profile in the field lines associated with the return 
current. This means that both simulations have approximately 
the same magnetic shear profile along the PIL (neglecting the 
effect of different distant direct currents). From Equation (l24l> . 
it follows that both are expected to have very similar net cur¬ 
rent. This is what happens as inferred in Figure [7] from the 
curves of direct current, return current, and neutralization ra¬ 
tio of SP2 which match those of SP1. 

Finally, we emphasize that the above relationship between 
net currents and magnetically sheared PIL is significantly 
different from the Lor entz-force-driven shear discussed by 
IGeorgoulis et all(1201 2l) . Indeed, in all our zero-/!MHD simu¬ 
lations, there is no dynamical compression that could generate 
a Lorentz force which would shear the PIL. On the contrary, 
the magnetic shear along the PIL is caused by the motions im¬ 
posed in the photosphere, and this magnetic shear generates a 
force-free net current. 

5.3. Net currents versus neutralization predictions 

The possibility of generating coronal magnetic fiel ds carry- 
ing_a net current is at vari ance with the conjectures of lMelrose i 
(fT99TT) and lParkerl d 19961) who argued that twisted and sheared 
coronal fields should be perfectly current-neutralized. The 
main difference with our results relies in the fact that both au¬ 
thors built a conjecture using 2.5D considerations which were 
then directly transposed to 3D geometry. The cornerstone of 
their conjectures is that any 3D twisted/sheared magnetic field 
embedded in a potential field or field-free environment can be 
equivalently described by a flux tube connecting two parallel 
planes and set in the same (but 2.5D) magnetic environment. 

The above assumption is valid for any 3D flux tube whose 
two photospheric magnetic polarities are fully separated ei¬ 
ther by an ambient potential field or a field-free environment, 
i.e., when there is no magn etic shear at the PIL. This is the 
implicit assumption of lMelrosei ( 1991 ) who limited his con¬ 
siderations to the case of twisting and shearing motions not 
extendin g to the PIL . The same assumption is also implicitly 
made by lParkerl d 1996 ) when using a cylindrical twisted flux 
tube of finite radius as a model for an element of an AR mag¬ 
netic field which is fragmented in the photosphere. In such 
cases, a 2.5D analysis shows that electric currents should be 
neutralized (se e Ap pendixlAt. This agrees with our 3D deriva¬ 
tion of Section l5/2l in which we demonstrated that full current 
neutralization occurs when a PIL is fully embedded in a po¬ 
tential field or field-free environment. 


On the contrary, when a 3D flux tube has magnetic shear 
along its PIL, the above cylindrical description of an AR is 
not valid. Indeed, an electric current is then flowing along the 
PIL where the two pola rities are in con tact. Then, t he cyl in- 
drical approximation of lMelrose ( 1991 ) and lParkerl ( 1996 ) is 
not relevant. As a consequence, the associated conclusion that 
currents should be neutralized is not applicable. A fully 3D 
analysis is then required to predict whether or not currents 
should be neutralized. Such a 3D analysis actually show that a 
net current should be expected when magnetic she ar is present 
along the PIL of the 3D flux tube (cf. Section [L2l i. 

6. CONCLUSIONS 

In this study, we used 3D MHD numerical simulations to 
analyze the properties of electric currents in line-tied coro¬ 
nal fields generated by photospheric flows in bipolar ARs. 
We showed that typical photospheric flows, such as twist¬ 
ing and shearing motions, invariably produce both direct and 
return currents in the line-tied coronal fields. We find that 
these photospheric flows can create both neutralized and non- 
neutralized currents. 

Using Ampere’s law, we provided a physical origin to the 
build-up of force-free net currents in coronal magnetic fields. 
They arise fro m t he development of magnetic shear along 
PILs (Section 15.2b . Th is c onclusion agrees with the inde¬ 
pendent study of ITorok et al] (120141) who showed that mag¬ 
netic flux emergence can also produce a net coronal current 
which simu ltaneously develops with magnetic shear at the 
PIL. ITorok etlikl ( 2014 ) and our study thus set a theoretical 
framework for understanding the properties of electric cur¬ 
rents in ARs. They both show that, in general, net currents 
can be formed in the corona from various, independent and/or 
combined processes: e.g., magnetic flux emergence, photo¬ 
spheric flows, and by extension, any mechanism that can gen¬ 
erate magnetic shear along a PIL. 

Net currents in ARs can therefore develop in a large 
variety of cases. On the contrary, the production of perfectly 
current-neutralized magnetic fields in the solar atmosphere 
is a special case that requires a rather uncommon condition: 
a main PIL without magnetic shear. This is unlikely for 
both emerging and evolved/decaying ARs. Indeed, sev¬ 
eral observational and numerical studies show that newly 
emerged ARs generally possess a strongly s heared PIL (e.g., 
Manchester et al] l2004t ICanou et al.1 l2009t IGeorgoulis et all 

2012t iToriumi et all 120131: iPoisson et al.l 120151) . Next, 


evolved/decaying ARs are often associated with the presence 
of Her filaments, which are cold dense material supported in 
highly shear ed/twiste d magn etic fields lying_above PI Ls (e.g.. 


van Ballegooiien & Martens 1989; Aulanier & Demouiin 

1998 

: Gibson et al. 2004 Schmieder et al.l ;2006 ! ; Jing et al. 

201C 

; Mackay et akl[20l(I Xia et al. 201jjT 


Furthermore, starting from a configuration with an un¬ 
sheared PIL, it will remain so if the shear component of pho¬ 
tospheric flows is not extending to the PIL. This is also un¬ 
likely since the opposi t e is typically obse r ved in AR s (e.g. , 
Vemareddv et all 1201 2t IWang et akl 120121: iGuo et al] 120131 : 
Liu et al.l 120131) . In fact, several observations suggest shear¬ 
ing profiles that would be analogous to our strongly non- 
neutralized run S4 (as inferre d from the pho t ospheric motion s 
of magnetic polarities; e.g., ISu et alJl2007t ISun et al.ll20T2 ). 
Finally, the absence of significant shear is more likely to hold 
for evolved ARs possessing isolated magnetic polarities (i.e., 
magnetic polarities far away from each other). 

The two main sources of AR currents, emergence and hor- 
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izontal photospheric motions, would thus be expected to pri¬ 
marily produce net coronal currents. Nevertheless, a few ob¬ 
servational studies seem to indicate that ARs with neutralized 
currents may be as numerous as ARs w ith a net current ( e.g., 
iWilkinson et al.ll992tlWheatlandi2000i) . If we consider a typ¬ 
ical magnetic field of R max = 2000 G for newly emerged ARs, 
and 15 Mm for one spatial unit of our non-dimensionalized 
simulations, the strength of the net current in our strongly 
non-neutralized magnetic fields can reach ~ 0.7-4 x 10 12 A. 
If the magnetic field is scaled to B max = 200 G for de¬ 
caying ARs, the strength of the net current ranges between 
~ 0.7 - 4 x 10 11 A. Since the actual measurement precision 
is about 10 11 A, it then remains to be proved that current- 
neutralized ARs are truly current-neutralized. This requires a 
systematic analysis of both the current-neutralization and the 
magnet ic shear at the PIL for a ny st udied AR, which w as not 
done in IWilkinson et alJ (119921) and lWheatlandl (120001) . Such 
dedicated studies could then be used to further test the con¬ 
clusions derived in the present paper. 

Finally, even though our MHD simulations systematically 
report the presence of return currents, net currents in coro¬ 
nal magnetic fields do exist. Therefore, eruption mod¬ 
els based on magnetic configurations possessing net cur- 


rents ( e.g 

. van Tend & Ku]3erus!|1978 

; Hevvt 

ierts et al.II 1982h 

Lin et al. 

19981 Titov & Demoulin 

1991 

Kliem & Torok 

2006; Demoulin & Aulanier 20lCj) are a simplified, but valid. 


description of pre-eruptive magnetic fields in ARs. Our re¬ 
sults that relate net coronal currents to magnetically-sheared 
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Fig. 11.— PDF of the fractional flux, f (Equation <B2> ). of the magnetic field (left) and of the electric current density (right), for the twisting simulation T4 at 
/ = 60 tA- Top: linear scale. Bottom: log-scale. 
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APPENDIX 


A. ELECTRIC CURRENTS IN CYLINDRICAL GEOMETRY 

Let us consider a twisted magnetic flux tube in cylindrical geometry. The cylindrical coordinates r, 9, and z, respectively 
describe the distance to the axis of the flux tube, the angle around its axis, and the position along its axis. Using the integrated 
form of Ampere’s law, the total current, /, flowing across a disk of radius r, is 

Bo f Kr) = 2 nr Bg(r ), (Al) 


where / (equals +1 or -1) is included so that I(r ) is a positive function in the flux rope core where the direct current is located. In 
the flux rope core, the current I is therefore a growing function of r. The region of return current is located where I is a decreasing 
function of r, so for 


d(rB g ) 

dr 


< 0 . 


(A2) 


The existence of return currents is thus simply constrained by the existence of r / > 0, such that 

C 

Bg(r ) < — , for r> n, 
r 


(A3) 


where C is an integration constant. The region where Bg decreases faster than 1/r is the region of return current. 
The flux rope current is fully neutralized if there is a finite radius, r c , such that I(r > r c ) = 0, i.e., 


Bg(r) = 0, for r > r c . (A4) 

From Equation ( 1A 1 b . it is straightforward to show that if a twisted magnetic flux tube is confined (i.e., if it has a finite radius 
R = r c ), then the total electric current carried by the flux tube vanishes. 

These conditions are valid regardless of the force-freeness of the magnetic field. 


B. SOLENOIDAL CONDITIONS WITH OHM 

As mentioned in Section [2Tl the solenoidal condition for the magnetic field is not numerically imposed in the code. However, 
to show that it does not affect the evolution of the magnetic field in our simulations, we compute the fraction of non-conserved 
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TABLE 2 

Statistics of the pdf of fractional flux, fi . 


X 

< / (X) > 

median of L/j(X)| 

< L/j(X)| > 

<T(f(X )) 

V</(X)r> 

B 

-1 x 1(T 10 

2.3 x 1(T 7 

2.8 x 1(T 5 

1.7 x 1(L 4 

1.7 x 1(T 4 

j 

2.1 x 1CT 11 

3.3 x IQ- 8 

2.9 x 1(T 7 

2.5 x 1(T 6 

2.5 x 1(T 6 


Note. — B and j respectively are the magnetic field and the electric current 
density. <> and cr respectively are the mean and standard deviation. Results are 
presented for the simulation T4 at t = 60 . 


flux, or fractional flux, /, within each cell, i, of the mesh, such that 

(Bl) 


(B2) 

where |X, | and V • X, are re specti vely the n orm of X, and its divergence, within the mesh-cell, i, of volume, v,, bounded by the 
surface, s, (e.g., IWheatland et al.ll2000HValori et al .112013l) . 

We compute the set of /(B) values for the twisted flux tube simulation, T4. The corresponding probability density function 
(PDF) and its statistics are displayed in FigureQT]and Table[2] As one can see, the PDF is well centered on zero. The two bumps 
in the wings of the PDF are contributions from the boundaries. 

Both Figure|TT]and Table[2]show a mean, a median, and a standard deviation of f, that are all smaller than ~ 10 4 . This means 
that the amount of non-conserved magnetic flux within each mesh-cell is typically ~ 10 4 times weaker than the local magnetic 
flux. Therefore, even though the solenoidal condition is not numerically treated within the OHM code, the non-conservation of 
magnetic flux remains very weak. These results are representative of the set of numerical simulations performed and presented in 
this paper. We thus conclude that the solenoidal condition is well enough verified to ensure that the very weak non-conservation 
of magnetic flux does not affect the evolution of the system in each of our line-tied simulations. 

We further compute the fractional flux for the electric current density (/(j)) to check that its divergence is indeed zero (as one 
would expect from applying the divergence operator to Ampere’s law). The associated PDF, and the statistics of the PDF, are 
displayed in FigureQT]and Table[2] As for the magnetic field, V ■ j = 0 is well preserved in our line-tied simulations. 

C. NEUTRALIZATION RATIO: TRANSITION PHASE 

The neutralization ratio curves of the partially current-neutralized cases all present a transition phase between t = 0 and 
t ~ 2Q1 a- This transition phas e is distinguished by two specific periods, t = [0; 10]f^ and t — [10; 20] 

As mentioned Section 13.21 the direct/return currents are computed by extracting the negative/positive current density at the 
photospheric positive polarity. During the early evolution of the system (when the driving is zero and/or extremely weak), the 
computed total direct/return current, and hence neutralization ratio, are all dominated by the noise. This noise is essentially 
due to the presence of magneto-acoustic waves at early times. While the initial potential field is analytically potential and at 
equilibrium, it is not numerically because of the discretization. This well-known effect leads to the generation of magneto-acoustic 
waves inducing compression of the magnetic field. Such a compression creates transitory neutralized currents in each magnetic 
polarity of the system because they are not associated with magnetic shear build-up at the PIL. This is why the neutralization 
ratio is initially well-defined and equal to 1 for all simulations. A precise analysis shows us that the currents generated by the 
photospheric motions progressively become dominant (i.e., their strength is ~ 100 times larger than the noise) for t > 1 Or,}. 
The transition from a system dominated by noise currents towards photospherically-generated currents is thus responsible for the 
early evolution of the neutralization ratio, up to 10f,i. 

The second transition appears between t = [10; 20] Lt, which corresponds to the main acceleration period of the temporal 
ramp function. This transition follows the evolution of the temporal ramp function. Such a transition could be produced by 
two competitive mechanisms: (1) non-force-free effects due to the f ast a ccel erati on of the photospheric velocities, and (2) the 
saturation of currents due to field line length (as discussed in Sections [3.2l and l4~2l i. 

D. RETURN CURRENT DECREASE IN SHEARING RUNS 

Let us consider a magnetic field line at a time t. We define p as the distance between its two photospheric footpoints. a'o and vo 
are the distance between both photospheric footpoints in the x and y direction respectively, such that 

xo=psin</> (Dl) 

yo=pcos(p, (D2) 

where <f> is the angle between the field line footpoints and the normal to the PIL at the photosphere. The variation of these 
distances during an infinitesimal time, dr, is 


f V-X; dVi 

Z(X) = LL- 

l IX,] d Si 
_ Vj V • X; 

~ Si |X, | ’ 


dxo = dp sin </> + p cos </> d</> 
dyo = dp cos </> - p sin <f> d<p . 


(D3) 

(D4) 
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The photospheric motions being solely applied in the x direction, it follows that dvo = 0. Combining with Equations (ID3b and 
dD4b . one obtains the following equations 

dp = sin 0 dxo (D5) 

pd0 = cos 0 dxo ■ (D6) 

Then, we define |e| = \4> - 0o| 1 • Replacing in Equations (ID5I) and (ID6I) and expanding terms to a first order in e, one obtains 

dp «(sin 0o + e cos 0o) dxo (D7) 

pd0 *(cos 0o - e sin 0o) dxo • (D8) 

When the footpoints segment is close to the normal to the PIL, 0o = 0. It follows that dp ~ edxo pd0 w dxo. The boundary 
driving essentially induces a rotation of the field line with regard to the vertical direction. In other words, it shears the magnetic 
field line, thus increasing its current density. On the contrary, when the footpoints segment is almost aligned with the PIL, 
0o = 7r/2. Then, dp ~ dxo » pd0 ~ -edxo. The boundary driving essentially increases the distance between the field line 
footpoints. Since this distance is related to the field line length (larger footpoints distance implies larger field line), it follows that 
the photospher ic motio n s ess entially increase the length of the field line, hence reducing its current density (cf Equation (IT6l >: 
see also lAulanier et al 120051) . 

We then conclude that the continuous shearing of each magnetic field line can generate two evolutionary phases for their 
current density: (1) a first phase of increase with magnetic shear, and (2) a phase of decrease due to a fast elongation of the field 
lines. The shortest magnetic field lines ( i.e ., with the smallest vo) should be the first to be affected by the decrease of current 
due to their elongation because they align more rapidly with the PIL. For a given xo, smaller values of vo induce larger values of 

tan0(f) = xo/yo, and hence, values of 0(f) closer to n/2 {i.e., the value for which the field line is aligned with the PIL). 




